Magnetotransport in disordered delta-doped heterostructures 
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We discuss theoretically how electrons confined to two dimensions in a delta-doped heterostructure 
can arrange themselves in a droplet-like spatial distribution due to disorder and screening effects 
when their density is low. We apply this droplet picture to magnetotransport and derive the expected 
dependence on electron density of several quantities relevant to this transport, in the regimes of weak 
and moderate magnetic fields. We find good qualitative and quantitative agreement between our 
calculations and recent experiments on delta-doped heterostructures. 
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I. INTRODUCTION 

Delta doping of semiconductor heterostructures has a 
history of leading to new and interesting physical phe- 
nomena, such as the fractional quantum Hall Effect fi 
This new physics has come from increasing the mobility 
of doped electrons by minimizing the disorder they feel 
due to the spatial separation between the dopants and the 
electrons. While spatial separation reduces the magni- 
tude of fluctuations in the potential due to dopant atoms, 
it is well known that disorder due to the random positions 
of donors in the delta-doped layer has important conse- 
quences for the transport properties of heterostructures. 
Recent efforts have attempted to gain a clearer experi- 
mental picture of the nano-scale distribution of electronic 
states in two dimensional electron gases (2DEGs)i^i^i^ 
Such electronic inhomogeneity on the nano-scale is be- 
lieved to be important for transport in a wide variety of 
strongly correlated electronic materials i^i^ 

A step towards a more systematic understanding of 
the effect of disorder on transport properties has been 
taken in recent experiments on delta-doped devices by 
Baenninger et aZfSilS where a series of small devices was 
fabricated in which the distance between the dopant layer 
and the 2DEG was varied in a controlled way. Some of 
these devices show very interesting resistance effects in 
a magnetic field. In particular, the dependence of the 
magnetoresistance on electron density has been suggested 
as evidence for charge density wave formationpiS 

In this paper, we argue that the observations of Ghosh 
et aZ i^^'^^d^d'^ are a manifestation of charge droplet for- 
mation in a 2DEG. To justify this point of view, we ana- 
lyze the charge distribution in a disordered 2DEG due to 
a delta doped layer, to demonstrate that in the devices of 
interest, the charge distribution likely consists of droplets 
of charge that sit in minima of the screened disorder po- 
tential. It is well known theoretically that the charge 
distribution in delta-doped heterostructures should have 
a droplet- like structure at low electron density.-^:-^ There 



has also been recent interest in electron droplets in the 
vicinity of the 2D metal insulator transition ; ^^d^ but rel- 
atively little attention has been paid to transport in in- 
sulating samples in a magnetic field (with the exception 
of Efros and co-worker a^^d^i^O ) . 

We derive expressions for the physical parameters of 
electron droplets in the non-linear screening regime that 
is relevant to the experiments of interest here, and then 
apply this picture to describe the transport in these de- 
vices. We observe that our picture implies that the tun- 
neling between droplets decreases as a function of mag- 
netic field in a manner that is consistent with experiment. 
The picture for the magnetoresistance at small fields is 
along the lines of that proposed by Glazman and Raikh;^ 
while at larger fields, the magnetoresistance crosses over 
to the behavior expected by Shklovskii and Efrosi^i^ 
These approachea^LSSii^ provide the magnetic field de- 
pendence of the resistance, but the picture of electron 
droplets yields non-trivial predictions for the dependence 
of resistance on parameters other than magnetic field, 
such as electron and dopant densities. In particular we 
show that in the regime of magnetic fields where the resis- 
tivity p varies with magnetic field B as p{B) cx exp[ai?^], 

_ 3 

that a c>c Tie ^ , where Ue is the electron density in the 
2DEG, even though the average tunneling distance be- 
tween droplets is much larger than the average inter- 
electron spacing. We also expect a to be temperature- 
independent at low temperatures, as observed in experi- 
ment. 

This paper is structured as follows: in Sec. Hi] we in- 
troduce our model for the delta doped heterostructure 
and elucidate the mechanism by which electron droplets 
form, then characterize how their properties depend on 
sample parameters. In Sec. IIIII we discuss the expecta- 
tions for magnetotransport, based on the droplet picture, 
and compare the implications with actual transport mea- 
surements. Finally, in Sec. IIVI we conclude and discuss 
the implications of our results for experiment. 
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Figure 1: Bandstructure of the heterostructure indicated by 
bold solid lines. The 2DEG is formed at the interface 2 = 
of AlGaAs and GaAs. The dashed lines indicate the random 
positions of the bottom of the conduction band in the 2DEG; 
the randomness arises from a spatially nonuniform distribu- 
tion of 5— dopants. The potential at the interface is due 
to the gate voltage Vg, the ionized dopants' potential and 
screening by conduction electrons in the interface. Ec and 
are the conduction and valence band energies, respectively, in 
bulk GaAs and Ep is the Fermi energy. 

II. MODEL 

In this section we introduce the physical situation that 
we are interested in, and the model we have of this sys- 
tem. We then derive in detail how disorder and screening 
effects give rise to a droplet-like electron density distri- 
bution in the 2DEG, and give the physical parameters of 
the electron droplets. We finally present a brief summary 
of the electron droplet picture. 

A. Main parameters 

The physical situation and several of the relevant pa- 
rameters for the physics we will consider in more detail 
later are summarized in Fig. [TJ 

There are two classes of parameters that control the 
physics in the devices of interest: those that are intrinsic 
to GaAs, and those that can be varied experimentally, 
through fabrication of a device, or use of a gate. The 
properties that are intrinsic to GaAs are k = 12.9, the 
dielectric constant of the GaAs or AlGaAs dielectric, and 
m — 0.067me, the effective electron mass in GaAs, which 
imply a Bohr radius for electrons in the 2DEG of as = 
AiTKeoh? /me^ ~ 10 nm. The following properties can be 
tuned in experiments: 

1. d, the distance between the gate electrode and the 
2DEG. Experimentally, d « 300 nm. 

2. A, the distance between the J— doping layer from 
the 2DEG. A = 10 - 80 nm. The typical A used in 
our calculations is 50 nm. 



3. Ud ^ 10^^ cm ^, the density of dopants in the S 
layer. 

4. Tie ~ 10"'^"'^ cm~^, the average density of electrons in 
the 2DEG in a typical sample. This will be the typ- 
ical value in our calculations. Some samples have 
lower electron densities, Ue ^ 10^°cm~^. 

5. Vg, the voltage at the gate electrode, which is 
metallic. 

In our theoretical treatment, we always assume that d 3> 
A, and Ud ^ n^- The dielectric constants for GaAs and 
AlGaAs are also assumed to be the same without loss of 
generality. 

B. Electron droplet formation in a 2DEG 

Our discussion follows similar lines to the work of 
Gergel' and Suris,^** who considered the formation of elec- 
tron droplets in a 2DEG when the average position of 
the bottom of the conduction band (see Fig. [T]) is above 
the Fermi energy Ep- However, our results differ quali- 
tatively from Ref. [3, and these differences are essential 
for connecting our results with experiment. Hence we 
present the effect of various experimentally relevant pa- 
rameters on the properties of the droplets in some detail. 

We assume that the donor positions are uncorrelated, 
and distributed with a white-noise Gaussian distribu- 
tion, and that the donors are ionized after compensat- 
ing the band-bending field, which provides the electrons 
for the two dimensional electron gas. We ignore effects 
due to incomplete and weak ionizationi^ For a mean S- 
doping density of rid = 10^^ cm~^ and a 2DEG-dopant 
layer separation of A = 50 nm, a complete ionization 
of the donors corresponds to a compensating potential 
Komp = eAnrf/(Keo) IV. Random fluctuations in the 
dopant potential cause the bottom of the conduction 
band in some regions of the 2DEG to lie below Ep; this 
leads to the appearance of electron droplets. In the rest of 
our discussion, we consider the potential (f) in the 2DEG 
after subtracting the compensating potential. Charge 
fluctuations about the mean can be of either sign. We 
also assume that the chemical potential is uniform across 
the sample. 

The charge density p{z,r) in the AlGaAs region is 

p{z,r) ^ en{r)S{z - X), (1) 

where n(r) = n_|_(r) — 7i_(r) is the surface charge density 
in the (5-layer and r is a two-dimensional vector in the 
(5-doping plane. The subscripts ± denote the sign of the 
charge fluctuation. The spatial distribution of the dopant 
atoms satisfies 

(n(r)n(r'))-(n)' = n,<5(r-r'), (2) 

and Ud = (n-|-)-|-(n_) is the total surface charge density of 
both polarities in the (5— layer. The dopant atoms create 
a fluctuating potential (/)(r) in the 2DEG, 
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(r) = Vcomp + 

1 



P -r , dr' / dz {p{z, r') - en^ [(j){r')] S{z)} ' 

ineoK J Jo 



(3) 



^(r-r')2 + z2 ^(r-r')2 + (2d~z)2 



In Eq. dS]), ne[(j){r')] is the electron charge density in the 
2DEG and the first and second terms under the square 
root signs represent direct and image contributions from 
the charge distribution respectively. The position of the 
image charges is to the left of the metallic gate electrode 
z — d in Fig. [TJ 

When the gate voltage Vg is small, the number of elec- 
trons rie at the interface will be small. Let us first esti- 
mate the mean square value of the potential fluctuations, 
{Scfy^}, at the interface by ignoring screening effects of fi- 
nite electron density, rig. Using Eqs. ^ and ^ one has 



(47rKeo)^ 



dr 1 



1 



Vr2 + A2 



^r2 + (2d-A)2 
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A(2d- A) 



(4) 



When the gate voltage Vg is increased, the bottom of the 
conduction band approaches the Fermi energy Ep- If po- 
tential fluctuations due to the dopants are such that in 
some regions of the interface, the bottom of the conduc- 
tion band now lies below Ep, a finite ne(r) will be found 
at such places (see Fig. [T]). 

Next we estimate the effect of screening on the mean 
square fluctuation ((J^^) due to a local electron density 
Ueir). Consider a region of size R in the (5-layer. From 
Eq. ^ , one can estimate the magnitude of fluctuations 
of n(r) in this region: 



7ri?2 



dr'((n(r)n(r')) - (nf) = 



7ri?2- 



(5) 



Thus ^/ {5'n?{R)) ~ nj jix^l'^R is the characteristic fluc- 
tuation of the dopant charge densi ty. If the local electron 
charge density, ?T.e(r), exceeds \J {5rP-(T()) , then poten- 
tial fluctuations due to these fluctuations in charge den- 
sity will be screened by a redistribution of the electronic 
charge Uf,. If the l ocal electron charge density is much less 
than -y/ ((5n2(i?)) , then potential fluctuations in a region 
of size R are not screened. Therefore, for a given electron 
density ng, we only need to take into account fluctuations 



in regions of size R < R, 
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limit of integration over r in Eq. ^ can be replaced with 
e] consequently 
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A(2rf- A) 

A2+i?2 



/(2d- A)2 + i?2x 1/4' 



(2d-A)2 + i?2y ■ y A2 + i?2 



(6) 



Equation ^ yields two simple limiting regimes for the 
relevance of screening effects: 



Ude^ 



iln (1 



i?c > 2d > A 
2d > Rc, A 



(7) 



The physical picture behind Eq. ([7]) is as follows. The 
metal electrode screens potential fluctuations on a length 
scale of d. Thus, if the fluctuations of the dopant density 
give a screening radius Rc that is larger than c?, we may 
ignore the screening effects of a finite electron density Ug- 
For the devices we are interested in, the important limit is 
generally where d is considerably larger than A, although 
Rc may be of the same order as, or greater than d at very 
low electron densities. We also note that Rc cannot fall 
below , since the Gaussian approximation [Eq. ([2])] 

that we used to determine it will no longer be valid. 

Let us now estimate the spatial dimensions of the elec- 
tron droplets localized in the minima of the smoothly 
fluctuating 2DEG potential. We first do this for moder- 
ate electron densities, for which we can assume 2fi S> Rc, 
to illustrate our logic, and then also consider the situation 
2d ^ Rc, where one has to deal with Eq. © numerically. 

From Eq. ([7]), it is clear that the confining potential 
increases as the distance to the (5-layer, A, is decreased. 
Nevertheless, the kinetic energy cost of electron confine- 
ment in the z-direction means that the localization length 
in the z-direction does not vanish. We assume that the 
electron wavefunction is localized in the z-direction in 
the GaAs substrate with a characteristic lengthscale zq. 
Hence, the kinetic energy of the electron is of the order of 
n^h'^Tr'^ / {2mzQ), (n is a natural number) and the typical 
distance of the electron from the dopant layer is of the 
order of A -I- zq. The potential fluctuation expression in 
Eq. (O then needs to be modified: 
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A + zo 



(8) 



We now minimize the total energy with respect to zo,n to 
get the binding energy En for the n**^ sub-bandji^ where 



I'K^I'^nc- Thus the upper En 



,1/2^2 



2mZo „ 40rKeo 
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Rc 



A + zo,„ 



(9) 
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50 
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Table I: Values of the screening radius Rc, the bound state en- 
ergies El and E2, and the penetration distance for the lowest 
sub-band 20,1 in the GaAs layer for different values of electron 
density rie. Also shown are the droplet sizes Rp, the highest 
energy A of the electrons occupying a droplet, and the local- 
ization lengths ^ for inter-droplet tunneling. The lengths Rc, 
zo,i, Rp and ^ are all in units of nanometers. 



electron density of ndocai will be of the order of 
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\ + zo 



(11) 



In an electron droplet of size Rp, the kinetic energy of the 
electron occupying the highest state in the droplet will 
be of the same order as the potential energy fluctuation. 
If Rp is small compared to A -I- zq (we shall see shortly 
that this is the case), the local fluctuation, Eq. pT|) is a 
linear-i? confinement. For such a confinement, the virial 
theorem suggests that the kinetic energy should be half 
of the potential energy. For the lowest sub-band, n = 1, 
we have 



A bound state is always possible with such a poten- 
tial because while the kinetic energy decreases as I/zq 
the magnitude of the potential energy decreases only as 
l/zo,n- The minimization leads to the following transcen- 
dental equation for 2;o,„: 




Eq. pOj) gives good results when 2d ^ Rc- For values 
of the electron density where this condition is not well- 
satisfied, then in Eq. Q, the more accurate expression 
for potential fluctuations, Eq. ^ should be used. Table|T] 
shows the numerically calculated values of Ei, E2, zo,i 
and zo^2 using Eq. with A replaced by \ + zq, for the 
potential fluctuations. 

The 2DEG electrons fill the lowest sub-band first and 
the second sub-band does not begin to be populated un- 
til the most energetic electron in the lower level reaches 
i?2. The basic result of this analysis is that droplets can 
form due to electrons filling the minima of the screened 
disorder potential. 

We now consider the dimension of the droplets in 
the plane of the 2DEG. The following arguments are 
valid for both zg > A and zq < A. In a region of 
size R <^ Rc, the charge density fluctuation is of the 
order of ri^J^ /n^^'^R ^ rig. The electrons begin fill- 
ing these small but deep regions first in an attempt to 
screen the strongest charge fluctuations and the density 
f^e, local — n^J^ /tt^/'^R of clcctrons in these droplets will 
be much larger than the mean 2DEG electron density rig. 
However, since the kinetic energy rises with confinement, 
this limits the density from becoming too large, since 
to high a density would prevent the formation of bound 
states. With an electron density of ne,iocai, correlations of 
the potential beyond the length scale R will get screened. 
Using the arguments we employed for {S(j)'^) earlier, we 
can similarly say that the potential fluctuation with an 
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A + zo,i 



(12) 



It should be noted that we can only make an order of 
magnitude comparison, as obtaining the exact magnitude 
of the screened Coulomb fluctuations where the electrons 
are pooled is beyond the scope of our discussion. To es- 
timate the wavevector fcmax we note that the droplet is 
effectively two-dimensional since motion deep into the 
GaAs substrate is not possible. Equating the number of 
droplet eigenstates, {k-^^^RpY /2, with the charge fluctu- 
ation Nc in the droplet, — ■n^^'^nj'^ Rp, (the number of 
electrons screening the potential fluctuation on a length 
scale R), we have fc^^x — 2n'J'^-K^I'^ jRp. Using this esti- 
mate for fcmax in Eq- (|12p . we arrive at the condition 
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'r 



p 



In 1 



Rn 



A + zo,: 



(13) 



and when we specialize to the case A -I- zo,i ^ as, the 
droplet condition, Eq. gives 



R. 



'.p ^ — \J 2aB(A + zo,i). 



(14) 



as the size of the droplet in the lowest sub-band. We note 
that in Ref. Q the form of the potential energy used to es- 
timate Rp was for the A <C i? limit, however, the conclu- 
sion of the analysis was that Rp ~ A. This invalidates the 
original assumption and Eq. (fT^ should be used instead. 
Equation (|14p gives for rig = 5 x 10^°cm~^ and zq.i ~ 
43 nm, Rp^i — 43 nm. This compares well with the nu- 
merical estimate, Rp ~ 50 nm (see Table HI. 

The estimate for Rp can be further improved by not- 
ing that Zq itself varies with Rp. Physically, the confining 
potential is weaker at smaller values of Rp, therefore zg 
increases as Rp is decreased. This effect will lead to an 
enhancement of Rp. For this purpose, one needs to min- 
imize Eq. ^ for Zq as a function of Rp <C (A + zq), 
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and then use the i?p— dependent zq in Eq. (jl3p to obtain 
the droplet size. We have not followed this quantita- 
tively more accurate but tedious procedure. Many of our 
quantitative estimates will be affected by our inability, in 
the present paper, to get a better estimate for Rp. How- 
ever our qualitative predictions, and in particular, the 
dependence of quantities on experimental parameters is 
not affected. In our quantitative estimates we will, for 
the sake of greater accuracy, use the numerically com- 
puted values rather than the analytic expressions that 
are strictly valid for 2d^ Rc- 

The density of electrons in the droplets, rie, local, is a 
little greater than the average electron density, Ug. For 
example, for rie = 5 x lO^^cm"^, we have rte, local — 
n^l"^ I'K^I'^Rp ~ 1.1 X lO^^cm^^. The local variation in 
the electron density is due to the movement of charge in 
the potential landscape from the "hills" to the "lakes" . 
The number of electrons in a droplet, iVg = ■K^^'^n'J'^ R^^ 
is only weakly dependent on the average electron density 
through zg- Extra electrons are accommodated in the 
2DEG through increasing the density of droplets. The 
separation between the droplet centers. 



^2./RJi;, (15) 

decreases with increasing average electron density rif, . We 
note that a previous discussion of droplet formation in 
2DEGs'''^ states that kp — Rc for the purposes of calcu- 
lating the conductivity. We correct that statement here, 
as our analysis makes clear that this is not the case for 
the regimes of interest. 

Next we discuss the height of the barrier separating 
two droplets, ii^barrier- This barrier is not simply the 
size of the fluctuation from the bottom of the poten- 
tial well, e-\/ ((50^), which is what one would find for 
"empty" potential wells. Instead we need to consider 
the occupation of the wells, and hence -Ebanior is given 
by the difference between the magnitude of the bind- 
ing energy E„ for confinement in the z direction and 
A = h^kl^^j2m = h^TT^/'^nJ'^ /{mRp), the highest en- 
ergy for the electrons occupying a well, reckoned from 
the bottom of the well. The inter-droplet barrier height 
^barrier is hcucc given by i^barricr = |£^i| " A. Addi- 
tionally, the typical number of electrons in a droplet 
is = -K^/^nJ^Rp ~ 9 for rie = 5 X IQi" cm'^ and 
Rp = 50 nm. Thus the mean level spacing in the droplets, 
5 ~ A/A^e = h^/{mRl) ~ 5 K. As A is typically of the 
order of the sub-band spacing \Ei — the second sub- 
band in the z direction is also likely to be populated to a 
certain degree, although at low enough electron densities 
and temperatures, this should not be a major concern, 
and due to the uncertainty of our estimates, we shall as- 
sume that only the lowest sub-band is populated. We 
do not expect this to significantly alter our conclusions. 
The distance r between the surfaces of two neighboring 
droplets is 

(16) 




(b) 



(r) 




Figure 2: Schematic picture of droplets. (a) A typical 
droplet (shown in solid shading) has a radius of the or- 
der of the distance A between the 2DEG and the 5— layer. 
For a dopant density Ud, the total number of electrons in 



the droplet is A^e 



,1/2 



Rp. For an average 2DEG electron 



density ne, the "catchment area" (shown in line shading) 



,1/2 



Rp/rie. 



RcRv 



the order of yJRcRp. (b) Cross-section of energy profile be- 
tween droplets. The electrons fill up to an energy A from the 
bottom of the potential wells. T he magnitude of the potential 
fluctuations is \Ei \ ~ ey^{S(j)'^) > A. 



The localization length for inter-droplet tunneling can be 
obtained from the size of the barrier, 



e = VV2mi;barricr = h/ ^2m{\Ei\ ~ A) . 



(17) 



Table H] lists the calculated values of A and ^ for different 
electron densities. Note that the localization length here 
is of the order of the Bohr radius in GaAs. This should 
be regarded as a coincidence. 



C. Summary 



r = hp — 2Rp — 2(y\J RcRp — Rp). 



In the previous section (Sec. Ill Bp . we showed that for 
low electron density delta doped heterostructures that 
give rise to 2DEGs, it is very natural to expect that elec- 
trons will organize themselves into droplets of charge, 
centered on the minima of the screened disorder poten- 
tial, for a quite reasonable range of parameters. The 
important length and energy scales of the droplets that 
relate to the transport properties of the device are the 
inter-droplet spacing, and the energy barriers between 
droplets, as summarized in Fig. [21 
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III. MAGNETOTRANSPORT 

In this section we consider several aspects of magneto- 
transport. In particular, we first consider the regime of 
low magnetic fields, where, due to quantum interference 
effects, one generically expects negative magnetoresis- 
tance, and we then move to higher magnetic fields, where 
one expects positive magnetoresistance due to shrinking 
of the localization length with increasing magnetic field. 
Our results are in good agreement with experiment. 



A. Small fields: negative magnetoresistance 

At small magnetic fields, many 2DEGs show negative 
magnetoresistance, i.e. the resistance decreases with in- 
creasing magnetic field, which is generically due to the 
magnetic field suppressing quantum interference of dif- 
ferent electron paths. The samples in Ref. [l^ show 
quadratic negative magnetoresistance at small fields: 



[7^(B)-7^(o)]/7^(o) = -(B/B,f 



(18) 



with Be — 0.15 T. In the same magnetic field range, 
the temperature dependence of the resistance obeys an 
Arrhenius law for temperatures greater than about 1 K. 

There are two contexts in which negative magnetore- 
sistance in insulators has been studied in the literature. 
The first is hopping conductivity in dirty semiconductors, 
and the second is tunneling between droplets of charge, 
similarly to the picture outlined in the previous section. 

In dirty semiconductors, electrical conductivity oc- 
curs due to hopping between different impurity sites 
and negative magnetoresistance is a consequence of sup- 
pression of the destructive interference of the Aharonov- 
Bohm phases acquired different trajectories j^^'^^ *^^ Such 
a theory gives a negative magnetoresistance cx \B\, 
the perpendicular magnetic field, for fields such that 
{BD^^^e/^)/cj)o > e-"/^^ K is the localization 
length, 00 = h/e is the fiux quantum, and D is the 
hopping distance.] For smaller magnetic fields the nega- 
tive magnetoresistance is quadratic in _B. If we use the 
value of ^ and D = r from our droplets scenario, this cri- 
terion suggests the negative magnetoresistance will be 
quadratic in the applied field up to _B ~ 0.2 T for 
rie = 5 X 10^° cm~^. 

In our case, the physical situation is somewhat different 
because instead of hopping between impurity sites as in 
a dirty semiconductor, the electrons in the 2DEG hop 
between droplets. The distance between the centers of 
neighboring droplets, lip = 2yJ R^Rp, is comparable with 
the droplet diameter 2i?p, unlike a dirty semiconductor 
where the electrons are localized at point-like impurity 
sites. This situation was studied in Ref. [U where it was 
shown that the resistance between two droplets behaves 
as 



7^(o) 



= (S/So)' 



1 



ne(cm ^) 


Bo(T) 


Bi(T) 


_D(nm) 


a(T-^), 
D = (3«?,r/4)i/3 


a(T-^), 
D = r 


10" 


1.45 


0.22 


33 


0.48 


9.6 X 10"* 


5 X 10^° 


0.42 


0.12 


95 


5.7 


0.84 


2 X 10^° 


0.19 


0.05 


180 


29 


13 


10^0 


0.10 


0.024 


271 


94 


62 



Table 11: Magnetic fields Bq and Bi as defined in Eq. (|20|l 
for different electron densities. These fields determine the 
crossover between negative and positive magnetoresistance. 
Also shown are the inter-droplet tunneling distance D = 
(3Zfp/4)^/^ and a, both of which appear in magnetoresistance 
expressions m Sec. IXITBI D = (3Zipr/4)i/^ is the tunnel- 



ing distance obtained by reconciling Eqs. (|19|l and (|20|) with 
Eq. (|21|l . Values of a axe also shown for comparison by as- 
sumimg the tunneling distance D is equal to the distance r 
between the surfaces of neighboring droplets. 



where Bq and Bi have the following approximate expres- 
sions in terms of our droplet parameters: 



-Bo ~ (l)o/{TTyalip), 
B, ^ 20o/(7r/?J, 



(20) 



cosh\B/Bi) 



(19) 



and we note that lip — 2^jRcRp, and we estimate that 
the spread of the wavefunction under the barrier in the 
direction perpendicular to the tunneling is j/o ~ 
Both Bq and Bi are determined by the field for which a 
fiux quantum is enclosed within an area that is of signifi- 
cance in the droplet picture. Bq is the field that encloses 
a fiux quantum in an area of the order of the area en- 
closed by the electron wavefunction tunneling under the 
barrier. Bi is the field that encloses a flux quantum in 
an area of the order of the "catchment region" . [See Fig. 
[5] . Physically, Bi is the field below which interference ef- 
fects are significant. If Bq > Bi, then at small fields, the 
magnetoresistance will be negative as in Eq. psp . with 
= Bf^ - Table H] lists Bq and Bi for various 
electron densities. 

For Tie = 5 X 10^" cm"^ we can infer from Table HIl 
that Be — 0.12 T, which is close to the experimentally 
observed Be — 0.15 We note that recent experi- 

ments on a quantum dot lattice^^ also displayed negative 
magnetoresistance that appeared to be well explained by 
the Glazman and Raikh approach. 

The positive magnetoresistance in Eq. (|19p arises from 
the shrinking of the localization length in a strong mag- 
netic field. This has been extensively studied in the 
contest of dirty semiconductors. In dirty semiconduc- 
tors, the 6*^^/^°-' dependence of magnetoresistance is 
predicted to cross over to a e^^^^ dependence at a high 
enough field. This is seen in the experiments we are dis- 
cussing too, while Eq. (|19p makes no such prediction. 
The high field data are better described by Shklovskii's 
expressions (see Sec. lIII B|) for magnetoresistance for tun- 
neling in dirty semiconductors , xhis is the subject of 
the following section. 
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B. Larger fields: positive magnetoresistance 

The important length scale for magnetotransport in 
the regime of magnetic fields where interference effects 
are negligible, is the length scale of the "unit cell" of 
the droplet array. What we mean is the regime Up > 
2lo, where hp — 2yJ RcRp is the separation between the 
centers of neighboring droplets, and Iq = ^Jh/eB is the 
magnetic length. 

Now, it is well known from Refs. [13 and[23l that once 
the magnetic field is large enough, B > Bi,so that quan- 
tum interference effects are negligible in comparison with 
the positive magnetoresistance arising from the reduction 
of the localization length, the magnetoresistance takes 
the form 



in the field range 



7^(o) ^ 



(21) 



(22) 



In Eq. (f2T|) . with D the typical tunneling distance as 
before, a is given by 



a 



(23) 



If we take a to be determined from Eq. , then we can 



(24) 



identify a ^ Bq^ , and the tunneling distance 



D 



The upper limit for B for the validity of Eq. pTjl . 
h/{e^D), is about 0.8 T for Ue = 5x lO^^cm-^. The quan- 
titative estimate for a is very sensitive to r and hp which 
limits its reliability. Using the extrapolation of Eq. 
one should observe 



A 



Rl 



A 

3/2 ' 



(25) 



where A — 2nJ Rp ^e^ /{n^/^fi^). A similar calculation, 
taking Shklovskii's expression for a, with D = r, also 

3 

yields a ~ A' /ni , where A' and A differ by a factor of 
order unity but have the same dependence on n^, Rp, and 
^. In both scenarios, the tunneling distance D increases 

1/2 

as l/rie as the electron density decreases, indicating 
that the rig dependence of a is robust; additionally, a is 
seen to be independent of temperature. This is exactly 
what was observed in Ref. [13, although in that case the 

_ 3 

Ue ^ dependence of a was ascribed to a Wigner crystal 
(or a charge density wave) where the density of localized 
states is expected to be equal to the density of electrons. 



We note that for the experimental parameters in Ref. [1, 
there is quantitative agreement between our prediction 
for a and the observed behavior. In particular, consider- 
ing D as given by Eq. ([24|) . we have for Ue ~ 10 
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cm 



-2 



a ~ 1.5 X 10^^m~^T-^/ni, and if rie ~ 5 x 10i°cm-^ 

3 

a ~ 6.4 X 10^^ m~'^T~^/ne . These values are about an 
order of magnitude of the values quoted in Ref. i9|, how- 
ever, we note as in Sec. Ill Bi that one of the main imped- 
iments to greater quantitative comparison with experi- 
ment is an inability to get a better estimate for Rp, which 
in turn affects our estimate of ^. If one assumes that the 
tunneling distance is equal to the distance between the 
neighboring droplet surfaces, D = r, then one obtains 

for ne ~ 5 X IQ^^ cmr'^ that a ~ 9.4 x lO^^ m-^T^V^I > 
which is closer to experiment. 

In our droplets picture such a behavior of D arises due 
to the rig dependence of the typical droplet spacing, and 
the dependence of a in Eq. (|25p is one of our main 
results. 

The expression for the magnetoresistance, Eq. (|2ip . 
is valid when D <C and for higher magnetic fields 

where D ^ 'o/C' (o^ equivalently, B ^ -^j^) it crosses 
over to another regimej^i^^ 



7^(o) 



where 



Bo = 



(26) 



(27) 



In this regime of magnetic field, the result of Ref. 
presented in Eq. (|19p is no longer valid. Note also that 
the magnetoresistance expression, Eq. (|26|) . rather than 
Eq. (PT]) or Eq. describes the correct behavior as 

the electron density is decreased. Our estimate for B2 
is 0.07 T when Z? ~ 95 nm. In the experimentSjiS a 
crossover from a quadratic to a linear magnetic field de- 
pendence (in the exponent) has been observed, and it 
is in qualitative agreement with the expected crossover 
B ^ -^j^ ^ 0.8 T, for the parameters mentioned above. 
Again, we are able to predict the dependence of this field 
on Ue] at low electron density. 



B2 oc 



en 



d 



This implies that the resistance should increase more 
quickly as a function of magnetic field in samples with 
lower rig, which is observed in experiments^'' 

The picture we have discussed in this section, is one in 
which the transport at increasing magnetic fields is domi- 
nated by barriers between droplets. As we have seen, the 
barrier transparency drops exponentially as the magnetic 
field is increased. This suggests that for large enough 
magnetic fields, individual droplets will be relatively well 
isolated from each other, and that one can think of the 
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droplets as an irregular array of quantum dots. This has 
further implications for transport which we will discuss 
in future work;^ 



IV. DISCUSSION 

The picture of electron droplets we have discussed here 
provides an appropriate framework for describing the 
phenomenology of the recent experiments in Ref. llOl . As 
we outline in Sec. lIIIBl we are able to use the droplet pic- 
ture to make both qualitative and quantitative compar- 
isons with experiment. [We note that some of the quanti- 
tative agreement may be fortuitous, however the depen- 
dence on 2DEG parameters should not be.] Importantly, 
we are able to reproduce the dependence of the magne- 
toresistance on the electron density. This strong agree- 
ment with experiment, in conjunction with our analysis 
of the electronic environment in the 2DEG, strongly im- 
plies that a picture of electronic droplets is more suitable 
to describe these, and similar experiments than a charge 
density wave scenario. 

One point that we did not discuss quantitatively, was 
how the non-linear screening is affected by the presence 
of a magnetic field. There has been some discussions 
of screening 2D electrons in a disordered potential in a 
magnetic field /^i^^'^° but this has tended to focus on the 
regime in which disorder is not too strong. We regard 
this as a very interesting problem that merits further in- 
vestigation; however, we note that experiment appears 
to provide some of the solution (hence our neglect of the 
issue here). Measurements of localized states in the quan- 
tum Hall regime^ that support a dot-like picture at low 
densities, find that the local electronic compressibility is 
essentially independent of the magnetic field, supporting 
our assumptions here. 

In our analysis we have assumed that the dopants are 
completely ionized. The degree of ionization is deter- 
mined by the internal field at the GaAs-AlGaAs junction, 
the external gate voltage Vg, and the temperature Tq at 
which the electron distribution over the dopants is frozen. 
Suppose the conditions are such that the donors arc not 
completely ionized. In particular, consudcr that the gate 
voltage Vg is zero and that the electron distribution over 
the dopants is frozen at a nonequilibrium temperature 
To which is less than the RMS potential fluctuation in 
the i5— layer. This case has been studied, for example, in 
Ref. m, where the authors calculate the dopant atoms' 
correlator, 



D(r 



(n(r)n(r')) — (n) 



(28) 



using a path integral approach. The probability of a fluc- 
tuation c(r) = n(r) — (n) is proportional to exp(— <I>[c]), 
where 

$[c] ^ — [ drc^M + / drdr'c(r)G(r - r')c(r'); 

2^^<i J KbTq J 

(29) 



G(r — r') is the interaction energy of two electrons in the 
(5— layer: 



G(r-r') = 



47renK 



1 



1 



|r-r'| |(r-r')2+4(d-A)2| 



(30) 



The dopant atoms' correlator is then given by 



D{r 



L'cc(r)c(r')e-*["1 



Dee 



Its Fourier transform D{q) can be shown to be 



g + go[l-exp[-2(7(d- A)]]' 



(31) 



(32) 



where go = n^e^ j (2eQKkBTQ) is the reciprocal Debye ra- 
dius. In Ref. [23, To — lOOK; for this value we may as- 
sume that qo ^ q. In the absence of dopant correlation, 
we have D{q) — Ud as in Eq.(l2]). For q{d — A) ^ 1, the 
correlator D{q) approaches a constant value. 



D{q) « no 



e2(d- A)' 



(33) 



We estimate hq ~ 2.5 x lO^cm^^ = 2 x lO^'^ii^. Note 
that the correlator in Eq. ((33|) has the same form as the 
uncorrelated case, Eq.®, except that the dopant den- 
sity rid is replaced by a smaller effective dopant density 
uq. This implies that the potential fluctuations in the 
2D electron layer (see Eq.(eq:meansqpotlimits)) are re- 
duced, and the characteristic length scale beyond which 
the donor charges are uncorrelated is much larger. The 
Rc corresponding to the reduced density is smaller than 
the original value by a factor of 20. Thus a much smaller 
electron density in the 2D layer suffices to screen the po- 
tential fluctuations from correlated dopants. The puddle 
size Rp is relatively less affected by dopant correlation. 
Lower values of Rc will take the system closer to metallic- 
ity since the puddles will merge once Rc becomes smaller 
than Rp. In the experiments we have studied, the ob- 
servations are better explained by assuming a complete 
ionization of the donors. 

We also find it worthwhile to say a few words on the 
effect of sample width on the puddles scenario. This has 
been studied numerically, for example, in Ref. [isl . The 
conclusion is that a reduction of the sample width di- 
minishes the effectiveness of screening of the potential 
fluctuations due to the dopant atoms, and the resulting 
enhanced potential fluctuations make puddle formation 
easier. In fact, the inter-puddle separation will increase. 
For studying ballistic transport in quantum wires, one 
should must work in a regime where the potential fluctu- 
ations are lower. This can be made possible by working 
with, say, gate voltages Vg where the dopants are incom- 
pletely ionized. 

It is interesting to view our results in the broader con- 
text of other physical phenomena observed in systems 
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with nanoscale electronic disorder. This includes other 
effects observed in low density delta-doped 2DEGs, and 
in systems such as the underdoped cupratesi^ The pic- 
ture of droplets of charge that merge as the doping is in- 
creased may also have relevance to the 2D metal-insulator 
transition.^2i2ii22 It is interesting to note that experi- 
ments on the insulating side of this transition have noted 
small jumps in the chemical potential as a function of 
Uer^ and we speculate that these small jumps may be 
relevant to the small amplitude oscillatory behavior of 
a as a function of rig observed in Refs. ^ and [l^- We 
note that our discussions here also lead to a plausible ex- 
planation of recent observations of a zero-bias anomaly 
in insulating 2DEG samples at low magnetic fields. ^^'"^^ 
This zero-bias anomaly might be due to exchange inter- 
actions between electrons on two closely coupled droplets 
having a Kondo effect similar to that observed in double 
quantum dot systems In the limit of higher density and 
larger magnetic fields, our picture should evolve contin- 
uously into existing droplet-based transport theory for 
the quantum Hall regimei^^i^i It would be interesting to 



try to connect our work with this scenario. Finally, our 
analysis of strong disorder in delta-doped heterostruc- 
tures could perhaps also be extended to a recent inter- 
esting proposal for creating strongly correlated electron 
systems by modulation doping near a heterojunction of 
two Mott insulators 
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